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' Abstract 



In [S. Effler, F. Ruskey, A CAT algorithm for listing permutations with a given number of inver- 
?H ' sions, I.P.L., 86/2 (2003)] the authors give an algorithm, which appears to be CAT, for generating 

permutations with a given major index. In the present paper we give a new algorithm for generating 
a Gray code for subexcedant sequences. We show that this algorithm is CAT and derive it into a 
CAT generating algorithm for a Gray code for permutations with a given major index. 



Q ; 1 Introduction 

u. 

We present the first guaranteed constant average time generating algorithm for permuta- 
tions with a fixed index. First we give a co-lex order generating algorithm for bounded 
compositions. Changing its generating order and specializing it for particular classes of 
compositions we derive a generating algorithms for a Gray code for fixed weight subex- 
cedant sequences; and after some improvements we obtain an efficient version of this last 
^ I algorithm. The generated Gray code has the remarkable property that two consecutive 
OO ■ sequences differ in at most three adjacent positions and by a bounded amount in these 
positions. Finally applying a bijection introduced in [7] between subexcedant sequences 
^ . and permutations with a given index we derive the desired algorithm, where consecutive 
^ I generated permutations differ by at most three transpositions. 

O ■ Often, Gray code generating algorithms can be re-expressed simpler as algorithms with 

the same time complexity and generating the same class of objects, but in different {e.g. 
lexicographical) order. This is not the case in our construction: the Grayness of the 
generated subexcedant sequences is critical in the construction of the efficient algorithm 
generating permutations with a fixed index. 

A statistic on the set &n of length n permutations is an association of an element of 
N to each permutation in 6„. For vr G &n the major index, maj, is a statistic defined by 
(see, for example, Section 10.6]) 



MAJ TT 



l<i<n 
7ri>7ri+i 

Definition 1. For two integers n and k, an n- composition of k is an n-sequence c 
C1C2 . . . c„ of non-negative integers with Yll=i — ^- -^^^ n-sequence b = 6162 ■ ■ - bn, 
is said b-bounded ifO<Ci< bi, for alH, 1 < z < n. 



1 



In this context 6162 • • • &n is called bounding sequence and we will consider only bounding 
sequences with either 6j > or 6j = bi-i = . . . = 61 = for alH, 1 < i < n. Clearly, 
6j = is equivalent to fix Cj = 0. We denote by C{k,n) the set of all n-compositions of 
k, and by C^{k^n) the set of 6-bounded n-compositions of k\ and if 6j > A; for all then 
C^{k,n) = Clk,n). 

Definition 2. A subexcedant sequence c = c\C2 ... is an n-sequence with < q < i — 1, 
for all i] and XlILi '^i called the weight of c. 

We denote by S{k, n) the set of length n and weight k subexcedant sequences, and 
clearly n) = C^{k, n) with 6 = 012... (n-l). 

2 Generating fixed weight subexcedant sequences 

We give three generating algorithms, and the third one generates efficiently combinatorial 
objects in bijection with permutations having fixed index : 

• Gen_Colex generates the set C^{k,n) of bounded compositions in co-lex order (de- 
fined later). 

• Genl_Gray which is obtained from Gen_Colex by: 

— changing its generating order, and 

— restricting it to the bounding sequence 6 = 01...(n — 1). 

It produces a Gray code for the set S{k,n), and it can be seen as the definition of 
this Gray code. 

• Gen2_Gray is a an efficient version of Genl_Gray. 

Finally, in Section HJ regarding the subexcedant sequences in S{k, n) as McMahon 
permutation codes (defined in Section |3]), a constant average time generating algorithm 
for a Gray code for the set of permutations of length n with the major index equals k is 
obtained. 

2.1 Algorithm Gen_Colex 

This algorithm generates C^{k, n) in co-lex order, which is defined as: C1C2 . . . c„ precedes 
did2 . . .dn in co-lex order if c„c„_i . . . Ci precedes dndn-i . . .di in lexicographical order. 
Its worst case time complexity is 0{k) per composition. 

For a set of bounded compositions C''{k,n), an increasable position (with respect to 
C^{k,n)) in a sequence C1C2--.Cn ^ C^{k,n) is an index i such that: 

• ci = C2 = . . . Ci-i = 0, and 

• there is a composition did2 . . .dn & C^{k, n) with q < di and q+i = (ij+i, Ci+2 = (ii+2, 

c = d 
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(^) 00000- 00001- 00002- 00003^ 

00013- 0023- 00033- 00133- 00233- 01233 



(b) 00000- 00003 

00033- 00233- 01233 

Figure 1: The path from the root to the composition 1 2 3 3 G C°'^'^^'^{9, 5): (a) before deleting 
redundant nodes (in boldface); and (b) in the generating tree induced by the call of Geii_Colex(9, 5) 
where redundant nodes are avoided. 



For example, for C°^^^'^(3,5) the increasable positions are underlined in the following 
sequences: 00010 and 20 0. Indeed, the first two positions in 00010 are not in- 
creasable since there is no composition in C°^^^^(3, 5) with the suffix 010; and the third 
position in 2 is not increasable because 2 is the maximal value in this position. 
Clearly, if £ < r are two increasable positions in c, then each i, i < i < r, is still an 
increasable position in c (unless 6j = 0). 

Here is the sketch of the co-lex order generating procedure for C^{k,n): 

• initialize c by the length n sequence ... 0; 

• for each increasable position i in c, increase Cj by one and call recursively the gen- 
erating procedure if the obtained sequence c is not a composition in C^{k,n), and 
output it elsewhere. 

The complexity of the obtained algorithm is 0{k) per generated composition and so 
inefficient. Indeed, too many nodes in the generating tree induced by this algorithm have 
degree one. Algorithm Gen_Colex in Figure [2] avoids some of these nodes. We will identify 
a node in a generating tree by the corresponding value of the sequence c; and a redundant 
node in a generating tree induced by the previous sketched algorithm is a node with a 
unique successor and which differs in the same position from its ancestor and its successor. 
For example, in Figured] (a) redundant nodes are: 0001, 000 2, 0013, 0023 and 013 3. 
These nodes occur when, for a given suffix, the smallest value allowed in an increasable 
position in the current sequence c is not 1, and this position is necessarily £, the leftmost 
increasable one. Algorithm Gen_Colex avoids redundant nodes by setting q to its minimal 
value e = k — Yl^^i (and Yl]=i "^^^ ^e computed for each i, 1 < i < n, in a pre- 
processing step). For example, in Figure [1] (b) there are no redundant nodes. However, 
in the generating tree induced by Gen_Colex there still remain arbitrary length sequences 
of successive nodes with a unique successor; they are avoided in procedure Gen2_Gray. 

Algorithm Gen_Colex is given in Figure [2] where i is the leftmost increasable position in 
the current sequence c, and r the leftmost non-zero position in c, and thus the rightmost 
increasable position in c is r if Cr < K and r — 1 elsewhere (6162 • • • &n being the bounding 
sequence). The main call is Gen_Colex(A;,ra) and initially c is 00 ... 0. (As previously, 
in this algorithm the function k 1— min{s | X]j=i — ^} can be computed and stored in 
an array, in a pre-processing step.) 

The induced generating tree for the call Gen_Colex(4,5) is given in Figure [3] (a). 
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procedure Gen_Colex(A:,r) 
global n,c,b; 
if fc = 
then print c; 
else if c[r] = b[r] 

then r := r — I ; 
end if 

£:=mm{s\j:;=im>k}; 
for i := £ to r do 

±fi = l then e /c — X]j=i ^[i] ' 
else e := 1 ; 

end if 

c[i\ := c[i] + e; 
Gen_Colex (A: — e,i) ; 
c[i] := c[i] — e; 
end do 
end if 
end procedure. 

Figure 2: Algorithm Gen_Colex. 

2.2 Algorithm Genl_Gray 

This algorithm is defined in Figure [Hand is derived from Gen_Colex: the order of recursive 
calls is changed according to a direction (parameter dir), and it is specialized for bounding 
sequences 6 = 012 ... (n — 1), and so it produces subexcedant sequences. It has the same 
time complexity as Gen_Colex and we will show that it produces a Gray code. 

The call of Genl_Gray with dir = produces, in order, a recursive call with dir = 0, 
then r — i calls in the for statement with dir equals successively: 

• 0, 1, ... 0, 1, if r — £ is even, and 

• 1,0,...1,0,1, if r-£is odd. 

In any case, the value of dir corresponding to the last call is 1. 

The call of Genl_Gray with dir = 1 produces the same operations as previously but in 
reverse order, and in each recursive call the value of dir is replaced by 1 — dir. Thus, the 
call of Genl_Gray with dir = 1 produces, in order, r — i calls in the for statement with 
dir equals alternatively 0, 1,0,.. ., then a last call with dir = 1. See Figure |3] (b) for an 
example of generating tree induced by this procedure. 

Let S{k,n) be the ordered list for S{k,n) generated by the call Genl_Gray(A;,n,0), 
and it is easy to see that S{k,n) is suffix partitioned, that is, sequences with the same 
suffix are contiguous; and Theorem H] shows that S{k,n) is a Gray code. 

For a sequence c, a > 1 and dir G {0, 1} we denote by first ( A; ; dir; c) and last(fc; dir; c), 
the first and last subexcedant sequence produced by the call of Genl_Gray(/c, r, dir) if the 
current sequence is c, and r the position of the leftmost non-zero value in c. In particular, 
if c = ... 0, then first(/c; 0; c) is the first sequence in S{k, n), and last(/i;; 0; c) the last 
one. 
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00000 



01210 

00120 ; 




00210 



00010 



01210 

00030 : 



00000 




00201 — 



00111 : 

01201 



(a) 

Figure 3: (a): The tree induced by the call of Gen_Colex(4, 5) with b 
by Genl_Gray(4,5). Terminal nodes are in bold-face 



(b) 

01234, and (b): that induced 
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Remark 1. 

1. For a sequence c, the list produced by the call Genl_Gray(fc, r, 0) is the reverse of the 
list produced by the call Genl_Gray(A;, r, 1), and with the previous notations we have 

last(/c; dir; c) = first(A;; 1 — dir; c), 

for dir G {0, 1}. 

2. Since the bounding sequence is 6 = 1 ... (n—l) it follows that, for c = ... QCj+i . . . c, 
Ci 7^ 0, first (A;; 0; c) is 

• aia2 . . . ai_iCiCi+i . . . c„ if /c < Y^jlXU " 1) = , where aia2 . . . aj_i is the 
smallest sequence, in co-lex order, in S{k,i — 1), 

• aia2 ■ ■ ■ ciiCi+i . . .Cn a k > liniHluHl ^ where aia2 ... is the smallest sequence, in 
co-lex order, in S{k + Ci, i). 



procedure Genl_Gray (fc,r,dir) 










global n,c,b; 










if fc = 










then output c; 










else if c[r] = r — 1 










then r :— r — 1; 










end if 










£ := min{s 1 "^7^^ >k}; 










. . J. (f-l)(«-2) . 
c • — ixi 2 ' 










if dir — 










then c[l]:=c[l]+e; Genl_Gray(A: - e,£,0) ; c[^] := c[ 


^]-e; 








dir :— (r — tj mod 2 ; 










for i :— £ + 1 to r do 










c[i] :— c[i] + 1 ; Genl_Gray (A: ~ 1 ,i , dir) ; 


dir := 


{dir 


+ 1) 


mod 2; c[i] :— c[i] — 1; 


end do 










else dir := ; 










for i := r downto ^ + 1 do 










c[i] := c[i] + 1 ; Genl_Gray (A; — 1 , i , dir) ; 


dir := 


[dir 


+ 1) 


mod 2; c[i] := c[i] — 1; 


end do 










c[£]:^c[e]+e; Genl_Gray(/c - e,£,l) ; c[£] := c[ 










end if 










end if 










end procedure. 











Figure 4: Algorithm Genl_Gray, the Gray code counterpart of Gen_Colex speciahzed to subexcedant 
sequences. 



Now we introduce the notion of close sequences. Roughly speaking, two sequences are 
close if they differ in at most three adjacent positions and by a bounded amount in these 
positions. Definition [3] below defines formally this notion, and Theorem H] shows that 
consecutive subexcedant sequences generated by Genl_Gray are close. 

Let s = S1S2 . . . Sn and t = tit2 ■ ■ - tn be two subexcedant sequences of same weight 
which differ in at most three adjacent positions, and let p be the rightmost of them (notice 
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that necessarily J9 > 3). The difference between s and t is the 3-tuple 

(ai, a2, as) = (sp_2 ^ Sp_i — tp-i, Sp — tp). 

Since s and t have same weight it follows that ai + a2 + as = 0; and we denote by 
— (ai, a2, as) the tuple (— ai, — a2, — as). 

Definition 3. Two sequences s and t in S{k,n) are close if: 

• s and t differ in at most three adjacent positions, and 

• if (ai,a2,a3) is the difference between s and t, then 

(ai, a2, as) G {±(0, 1, -1), ±(0, 2, -2), ±(1, -2, 1), ±(1, -3, 2), ±(1, 1, -2), ±(1, 0, -1)}. 

Even if the second point of this definition sound somewhat arbitrary, it turns out that 
consecutive sequences generated by algorithm Genl_Gray are close under this definition, 
and our generating algorithm for permutations with a given index in Section H] is based 
on it. 

Example 1. The following sequences are close: 01201 and 00301; 01003 and 01021; 00201 
and 01011; 01132 and 01204; the positions where the sequences differ are underlined. 
Whereas the following sequences are not close: 0211 and 01030 (they differ in more than 
three positions); 01201 and 01030 (the difference 3-tuple is not a specified one). 

Remark 2. If s and t are two close subexcedant sequences in S{k,n), then there are at 
most two 'intermediate' subexcedant sequences s', s" in S{k,n) such that the differences 
between s and s', between s' and s", and s" and t are ±(1, —1,0). 

Example 2. Let s = 010111 and t = 002011 be two sequences in 5(4, 6). Then s 
and t are close since they difference is (1, —2, 1), and there is one 'intermediate' sequence 
s' = 001 1 1 1 in 5(4,6) with 

• the difference between s and s' is (1,-1,0), 

• the difference between s' and t is (—1, 1, 0). 

A consequence of Remark [Tl2 is: 

Remark 3. If s and t are close subexcedant sequences and m is an integer such that 
both u = first (m; 0; s) and v = first (m; 0; t) exist, then u and v are also close. 

Theorem 4. Two consecutive sequences in S{k, n) generated by the algorithm Genl_Gray 
are close. 

Proof. Let s and t be two consecutive sequences generated by the call of Genl_Gray (A; , 0) . 
Then there is a sequence c = ciC2 . . . c„ and a recursive call of Genl_Gray acting on c (re- 
ferred later as the root call for s and t) which produces, in the for statement, two calls 
so that s is the last sequence produced by the first of them and t the first produced by 
the second of them. 

By Remark [TJl it is enough to prove that s and t are close when their root call has 
direction 0. 
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Let ^ and i ^ be the leftmost and the rightmost incrcasable positions in c (and 
■s,o c\ — ci — . . . = Cr-i = 0, and possibly = 0); and i and i + 1 be the positions where 
c is modified by the root call in order to produce eventually s and t. Also we denote 
m — k — Yl^=i s-nd e — m — ^'^^'"""^ . 

We will give the shape of s and t according to the following four cases. 

1. i — i and r — £ is even, 

2. i = i and r — £ is odd, 

3. i ^ i and the call corresponding to i in the for statement of the root call has direction 

(and so that corresponding to i + 1 has direction 1), 

4. i £ and the call corresponding to i in the for statement of the root call has direction 

1 (and so that corresponding to i + 1 has direction 0). 

Case 1. 

s = last(m — e; 0; 00 . . . CQ+i . . . c„) 

= first (m — e; 1; 00 . . . eQ+i . . . c„) 

_ r first(m - e - (£ - 2); 0; 00 . . . (£ - 2)eQ+i . . . c„) 
~ \ firstfm - e - (£ - 2); 0; 00 . . . (£ - 3)(e + 1)q+i 

and 

t = first(m-l;0;00...(Q+i + l)...c„) 

= first(m — e; 0; 00 ... (e — l)(c£+i + 1) . . . c„) 
= first(m - e - (£ - 2); 0; 00 . . . (£ - 2)(e - 1)(q+i + 1) . . . c^). 

Case 2. In this case s is the same as in the previous case and 

t = first(m-l;l;00...0(c£+i + l)...c„) 

_ f first(m-2;0;00...0(Q+i + 2)...Cn) if q+i + 2 < ^ 

~ \ first(m-e;0;00...0(e-l)(Q+i + l)...c„) if q+i + 2 > £ 

r first(m-e- (£-2);0;00...0(£-2)(e-2)(Q+i + 2)...Cn) if ci+i + 2 < £ 
\ first(m - e - (^ - 2); 0; 00 . . . (£ - 2)(e - 1)(q+i + 1) . . . c„) if q+i + 2 > ^. 

Case 3. In this case Ci — and 

s = last(m — 1; 0; 00 . . .Olcj+i . . . c„) 
= last(m-2;l;00...02Q+i...c„) 
= first(m-2;0;00...02ci+i...Cn), 

and 

t = first(m - 1; 1;00. . .0(ci+i + 1) . . .c„) 

r first(m-2;0;00...0(ci+i + 2)...c„) if q+i + 2 < i 
\ firstfm - 2; 0; 00 . . . 1(^+1 + 1) . . . c„) if a+i + 2>i. 



if e = e-i 

. . . c„) if e < £ — 1, 
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■ • • — 0000401 - 0003401 - 0023401 - 0123401 

Figure 5: Four successive q-terminal nodes in the generating tree induced by the caU Genl_Gray(ll,7,0) 
which generates the hst 5(11, 7). 

Case 4. As previously, Q = and 

s = last(m — 1; 1; 00 . . . OlCj+i . . . Cn) 
= first(m — 1; 0; 00 . . . OIq+i . . . c„), 

and 

t = first(m - 1; 0; 00 . . . 00(ci+i + 1) . . . c„). 

Finally, by Remark [3] it follows that in each of the four cases s and t are close, and the 
statement holds. □ 

As a byproduct of the previous theorem and Remark [TJ 2 we have 

Remark 4. If s = S1S2 • • • s„ and t = tit2 ■ ■ - tn are two consecutive sequences generated 
by Genl_Gray and p is the rightmost position where they differ, then S1S2 . . . Sp-2 and 
tit2 ■ ■ ■tp_2 are the smallest, in co-lex order, sequences in S{x,p — 2) and S{y,p — 2), 
respectively, with x = Si + S2 + ■ ■ ■ + Sp_2 and y = ti + t2 + ■ ■ ■ + tp_2. Remark that 
S1S2 ■ ■ ■ Sp-2 = tit2 . . . tp-2, and so X = y, if s and t differ in two (adjacent) positions. 

2.3 Algorithm Gen2_Gray 

Since the generating tree induced by the call of Genl_Gray contains still arbitrary length 
branches of nodes of degree one it has a poor time complexity. Here we show how 
some of these nodes can be avoided in order to obtain the efficient generating algorithm 
Gen2_Colex presented in Figure [6l 

A quasi-terminal node [q-terminal node for short) in the tree induced by a generating 
algorithm is defined recursively as: a q-terminal node is either a terminal node (node with 
no successor) or a node with only one successor which in turn is a q-terminal node. The 
q-terminal nodes occur for the calls of Genl_Gray (fc, r, dir) when k = ^^''"'^^ . See Figure [5] 
for an example. 

The key improvement made by Gen2_Gray consists in its last parameter p, which gives 
the rightmost position where the current sequence differ from its previous one in the list 
S{k,n), and Gen2_Gray stops the recursive calls of more than three successive q-terminal 
calls. Thus, Gen2_Gray generates only suffixes of the form Cp_2Cp_iCp . . . c„; see Table [T]for 
an example. Since two consecutive sequences in the Gray code S{k,n) differ in at most 
three adjacent positions, these suffixes are enough to generate efficiently S{k,n), and to 
generate (in Section H]) a Gray code for the set of length n permutations having the major 
index equal to k. 

Now we explain how the parameter p propagates through recursive calls. A non ter- 
minal call of Gen2_Gray produces one or several calls. The first of them (corresponding 
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procedure Gen2_Gray(A:, r,dir,p,u) 




global n, c, b; 




if fc = or (p — r) > 3 and k = ^''^~^'> 




then output (p, li) ; 




else if c[r] = r — 1 




then r :— r — 1 ; 




end if 




^:=min{s| '(7'' >fc}; 




e :- k ; 




if dir = 




then c[£] ~ c[£] + e; Gen2_Gray (fc - .O.p.it) ; c[^] := c[£] - e; 




dir := (r — mod 2; 




for i :— £ + 1 to r do 




c[i] c[i] + 1 ; Gen2_Gray (fc — 1 ,i, dir ,j , fc — 1 + c[i]) ; dir := {dir + 1) 


mod 2; c[i] := c[i] — 1; 


end do 




else dir :— 0; 




for i :— r downto £ + 1 do 




if i = r then q :— p; v := u; else g=:i + l; u := c[i + 1] + fc ; end 


if 


c[i] := c[i] + 1 ; Gen2_Gray(fc — 1 ,i, dir , ; dir := (dir + 1) mod 2; 


c[i] := c[i] — 1 ; 


end do 




a £ = r then q :— p; v := u; else q:—£+l; v :— c[£ + 1] + k; end if 




c[£] := c[£] + e; Gen2_Gray (fc - e,^ , 1 , ; c[^] := c[^] - e; 




end if 




end if 




end procedure . 





Figure 6: Algorithm Gen2_Gray. 



to a left child in the generating tree) inherits the value of the parameter p from its parent 
call; in the other calls the value of this parameter is the rightmost position where the 
current sequence differs from its previous generated one; this value is i if dir = and 
i + 1 if dir = 1. So, each call keeps in the last parameter p the rightmost position where 
the current generated sequence differs from its previous one in the list S{k,n). Proce- 
dure Gen2_Gray prevents to produce more than three successive q-terminal calls. For 
convenience, initially p = 0. 

The last two parameters p and u of procedure Gen2_Gray and output by it are used 
by procedure Update_Perni in Section |4] in order to generates permutations with a given 
major index; u keeps the value of Ci + C2 + . . . + Cp, and for convenience, initially u = 0. 

Even we will not make use later we sketch below an algorithm for efficiently generating 
the list S{k, n): 

• initialize d by the first sequence in S{k,n), i.e, the the smallest sequence in S{k,n) 
in co-lex order, or equivalently, the largest one in lexicographical orders, and c by 
... 0, 

• run Gen2_Gray(fc, n, 0, 0, 0) and for each p output by it update d as: 2] := c[p— 2], 
d[p - 1] := c[p - 1], d[p\ := c[p]. 

Analyze of Geii2_Gray 

For a call of Gen2_Gray(/c,r,(izr necessarily k < ^^^^y^^, and if /c > and 
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sequence 


P 


permutation 


sequence 


P 


permutation 


012100 




214356 


010012 


5 


536124 


010300 


4 


324156 


001012 


3 


635124 


001300 


3 


423156 


000112 


4 


135624 


002200 


4 


413256 


000022 


5 


235614 


011200 


3 


314256 


000004 


6 


345612 


012010 


5 


215346 


000013 


6 


245613 


011110 


4 


315246 


000103 


5 


145623 


002110 


3 


513246 


010 3 


4 


645123 


310 


4 


123546 


010003 


3 


546123 


001210 


4 


523146 


010021 


6 


436125 


010210 


3 


325146 


001021 


3 


634125 


010030 


5 


435126 


000121 


4 


134625 


001030 


3 


534126 


31 


5 


234615 


000130 


4 


134526 


000211 


5 


124635 


000040 


5 


234516 


001111 


4 


624135 


000220 


5 


124536 


010111 


3 


426135 


001120 


4 


524136 


002011 


4 


614235 


010120 


3 


425136 


011011 


3 


416235 


002020 


4 


514236 


011101 


5 


316245 


011020 


3 


415236 


002101 


3 


613245 


011002 


6 


516234 


000301 


4 


123645 


002002 


3 


615234 


001201 


4 


623145 


000202 


4 


125634 


010201 


3 


326145 


001102 


4 


625134 


012001 


4 


216345 


010102 


3 


526134 









Table 1: The subexcedant sequences generated by the call of Genl_Gray (4, 6, 0) and their corresponding 

length 6 permutations with major index equals 4, permutations descent set is either {1,3} or {4}. The 
three leftmost entries (cp_2,Cp_i,Cp) updated by the call of Gen2_Gray (4, 6, 0, 0, 0) are underlined, where 
p is the rightmost position where a subexcedant sequence differ from its predecessor. 
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• k < lnJ^I—?)^ then this call produces at least two recursive calls, 

• ^^^^^^'^^'^^ < k < "^^^"^^ ^ then this call produces a unique recursive call (of the form 
Gen2_Gray(A;',r, ■,■,■) , with k' = k — ilzlHlz^^^ which in turn produce two calls, 

• k = then this call is q-terminal call. 

Sine the procedure Gen2_Gray stops after three successive q-terminal calls, with a slight 
modification of Ruskey and van Baronaigien's [4J 'C^T' principle (see also [5j) it follows 
that Gen2_Gray runs in constant amortized time. 



3 The McMahon code of a permutation 

Here we present the bijection ip : S{n) — (5„, introduced in [7], which have the following 
properties: 

• the image through ip of S{k, n) is the set of permutations in (3„ with major index k, 

• is a 'Gray code preserving bijection' (see Theorem [6]), 

• r is easily computed from cr and from the difference between s and t, the McMahon 
code of a and r, if s and t are close. 

In the next section we apply ip in order to construct a list for the permutations in 
with a major index equals k from the Gray code list S{k,n). 

Let permutations act on indices, i.e., for a = aia2 ■ ■ ■ Cn and t = ri r2 . . . r„ two 
permutations in (5„, a ■ r = anCTj.^ ... ar^. For a fixed integer n, let k and u be two 
integers, < k < u < n, and define [[u, k]] G &n as the permutation obtained after k right 
circular shifts of the length-w prefix of the identity in (3„. In two line notation 



II 7-u _ f 1 2 ■ ■ ■ k k + 1 ■ ■ ■ u u + 1 ■ ■ ■ n 

\u-k + l u-k + 2 ■■■ u 1 ■■■ u-k u + l ■■■ n 

For example, in 65 we have: [[3, 1]] = 3124 5, [[3,2]] = 2114 5 and [[5,3]] = 34512 
(the rotated elements are underlined). 

Let ip '■ S{n) — i- &n be the function defined by 

^Pitlt2 ...tn) = [[n, t„]] ■ [[n - 1, tn-l]] ■ . . . ■ [[z, U]] [[2, t2]] ■ [[1, h]] 

= I1[[M.]]. 

i=n 

Lemma 5 ([7]). 

1. The function ip defined above is a bijection. 

2. For every t = tit2 . . .t^ e S{n), we have maj IlLntt^' ^4] = Z^Li ^i- 
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The first point of tlie previous lemma says tliat every permutation vr G (3„ can be 
uniquely written as IlLnll^' some U's, and the subexcedant sequence tit2 . . .tn is 
called the McMahon code of tt. As a consequence of the second point of this lemma we 
have: 

Remark 5. The restriction of ip maps bijectively permutations in S{k, n) into permuta- 
tions in &n with major index equals k. 

Example 3. The permutation 7r = 521643G can be obtained from the identity by 
the following prefix rotations: 

123456it^ 456123 561243 125643 i^ 251643i^ 521643 ftl:^ 521643^ 
so 

vr = [[6, 3]] -[[5, 4]] -[[4, 2]] .[[3, 2]] -[[2,1]] .[[1,0]], 

and thus 

MAj7r = 3 + 4 + 2 + 2 + l + = 12. 

Theorem [6] below states that if two permutations have their McMahon code differing in 
two adjacent positions, and by 1 and —1 in these positions, then these permutations differ 
by the transposition of two entries. Before proving this theorem we need the following 
two propositions, where the transposition (m, v) denote the permutation vr (of convenient 
length) with ti{i) = i for all except 7r('u) = v and 7r(w) = u. 

Proposition 1. Let n,u and v be three integers, n>3,0<u<n — 2, 1 < v < n — 2, 
and a,T E &n defined by: 

• a = [[n, m]] ■ [[n — l,f]], and 

• T = [[n,u + l]]- [[n-l,v-l]]. 
Then 

T = a ■ {n,v) . 

Proof. First, remark that: 

• [[n,u + 1]], is a right circular shift of [[^, m]], and 

• [[n — 1, f — 1]] is a left circular shift of the first {n — 1) entries of [[n — 1, v]], 

and so (y{i) = T{i) for alH, 1 < i < ra, except for i = n and i = v. □ 
Example 4. For n = 7, u = 4 and v = 3 we have 

. a = [[n, u]] ■ [[n - 1, v]] = [[7, 4]] ■ [[6, 3]] = 7 1 2 4 5 6 3, 

• r = [[n,u + l]] ■ [[n-l,v-l]] = [[7,5]] ■ [[6,2]] = 7134562, 

• {n,v) = (7,3), 
and T = a ■ {n,v) . 



13 



Proposition 2. If tc & (3„ and {u,v) is a transposition in ©„, then 

7i~^ ■ {u,v) ■ n = {n^^ (u) , n^^ (v)) . 

Proof. Indeed, (vr^-*^ ■ {u,v) ■ 7r)(i) = i, for all i, except for i = 7c^^{u) and i = n^^^v). □ 

Theorem 6. Let a and r be two permutations in ©„, n > 3, and s = SiS2 • • • s„ and 
t = tit2 . . .tn their McMahon codes. If there isaf,2<f<n — 1 such that ti = Si for 
all i, except tf = sj — 1 and = + 1, then r and a differ by a transposition. More 
precisely, 

where 

1 1 

«= n [[^'^*]]= n tt^'^^]]' 

and u = f + 1, V = Sf. 
Proof. 

• ^ = ntntt^.^d]. and so r ■ = nf=„[[^, ^^B, and 

• ^ = nLntt^' ^i]]' and a ■ = U(=n[[^^ ^i]]- 
But, by Proposition m 

/ / 

i=n i=n 

or, equivalently 

T ■ = a • ■(/ + !, s/), 
and by Proposition |2l the results holds. □ 

The previous theorem says that a and r 'have a small difference' provided that their 
McMahon code, s and t, do so. Actually, we need that s and t are consecutive sequences 
in the list S{k, n) and they have a more particular shape (see Remark H]). In this context, 
permutations having minimal McMahon code play a particular role. 

It is routine to check the following proposition (see Figure [7] for an example). 

Proposition 3. Let n and k be two integers, < k < "^•"^'''^ ; a = aia2 . . . a„ be the 
smallest subexcedant sequence in co-lex order with 'Y^^=i Q-j = k, and a = an,k = V'('^) be 
the permutation in having its McMahon code a. Let j = max{z : Oj 7^ 0}, that is, a 
has the form 

012...(j-3)(j-2)aj00...0. 

Then 

{j-Oj-i if 1 <i < 3 - {oj + 1), 
2j -aj-t if J - {aj + 1) < z < j, (2) 
i if i > j. 
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Figure 7: The permutation a = 21654387 with the McMahon code a = 1 2 3 4 3 0, the the smahest, 
in co-lex order, subexcedant sequence in S'(13,8), see Proposition |3l 

Remark 6. The permutation a defined in Proposition in [3] is an involution, that is 
qT^ = a. 

Combining Proposition [3] and Remark [6], Theorem [6] becomes in particular 

Proposition 4. Let a, r, s andt be as in Theorem\^ In addition, let suppose that there 
isaj,0<j<f — 1, such that 

1. Si = ti = for j < i < f — I, and 

2. ifj>0, then 

— Sj = tj 7^ 0, and 

— Si = ti = i — 1 for I < i < j . 

Then 

T = a-{<P,{f + l),Msf)) 

with 

{j — Sj — i if I < i < j — {sj + 1), 
2j-sj-t if J -{sj + 1) <t<j, (3) 
i if i > j. 

Notice that, the conditions 1 and 2 in the previous proposition require that siS2 . . . s/_i = 
tit2 ■ ■ - tf-i be the smallest subexcedant sequence, in co-lex order, in S{f — 1) with fixed 
value for X]f=i^ — J2i=i ^i- Also, for point 2, necessarily j > 2. 

4 Generating permutations with a given major index 

Let a and r be two permutations with their McMahon code s = S1S2 • • • Sn and t = 
tit2 . . - tn belonging to S{k, n), and differing in positions / and / + 1 by 1 and —1 in these 
positions. 
Let 

• V = Sf — tf E { — 1, 1}, and 
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If S1S2 . . ■ Sf-i is the smallest sequence in S{x, f — 1), in co-lex order, then applying 
Proposition m it follows that the run of the procedure transpCf , /, x) defined in Figure [8] 
transforms a into r and s into t. 



procedure transpC 




) 




J := inin|z : — 




x}; 




if u = 1 








then a := a ■ {(, 




-1), 


fMf])); 


else (T := a ■ {(, 




-1), 


i>Mf\ + ^))' 


endif 








s[f] ■= s[f] " « 








s[f + l] ■.= s[f~ 


)-i] + 


v; 




end procedure . 









Figure 8: Algorithm transp, where (pj is defined in relation ([3]). 

Let now / be the leftmost position where two consecutive sequences s and t in the 
list S{k,n) differ, and a and r be the permutations having they McMahon code s and 
t. By Remarks [2] and m we have that, repeated calls of transp transform s into t, and a 
into r. This is true for each possible 3-tuples given in Definition 12] and corresponding to 
two consecutive subexcedant sequences in S{k,n), and algorithm Update_Perm in Figure 
[3 exhausts all these 3-tuples. 

For example, if s and t are the two sequences in Example [2] with they difference 
(1, —2, 1), / = 2 and x = 0, then the calls 

transp (1, /, x) ; 
transp(-l, f + l,x + s[f]); 

transform s into t and a into r. 

Algorithm Gen2_Gray provides p, the rightmost position where the current sequence 
c differs from the previous generated one, and u = X^iLi'^*- Algorithm Update_Perm 
uses /, the leftmost position where c differs from the previous generated sequence, and 
^ = Ei=i Ci. 

Now, we sketch the generating algorithm for the set of permutations in Gn having 
index k. 

• initialize s by the smallest, in co-lex order, sequence in S{k,n) and a by the permu- 
tation in &n having its McMahon code s, 

• run Gen2_Gray(fc, ra, 0, 0, 0) where output is replaced by Update_Perm(p, u) . 

The obtained list of permutations is the image of the Gray code S{k,n) through the 
bijection defined in relation ([1]); it consists of all permutations in (5„ with major index 
equal to k, and two consecutive permutations differ by at most three transpositions. See 
Tabled] for the list of permutations in <Sq and with major index 4. 
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procedure Update_Perm(p, u) 


X :— u — c[p] — 


c[p-l]; 


if p - 2 > 1 and s[p - 2] = c[p - 2] 


then f '■= p — ^ 


L; 


else f := p — X 


I; x-~x- c[/] ; 


endif 




(ai,a2) := {s[f] 


-c[/],,s[/ + l]-c[/ + l]); 


if f + 2> n then 03 := 0; else 03 := s[/ + 2] - c[/ + 2] ; endif 


if ai > then t; := 1; else t; —1; endif 


case (01,02,03) of 


±(1,-1,0) 


transpCw, /, x) ; 


±(2,-2,0) 


transpCw, /, x) ; transpCu, /, x) 


±(1,-2,1) 


transpCw, /, x) ; transpC— u, / + 1, 2: + ; 


±(1,-3,2) 


transpCw, f,x); trsinspC— w, f + 1, x + s[f]) ; transp(— ti, / + 1, a: + s(/]) ; 


±(1,1,-2) 


transpCw, / + 1, x + s[/]) ; trcinspd;, /, x) ; transpCw, / + 1, x + s[/]) ; 


(1,0,-1) 


transpCl, /, x) ; transpCl, f +l,x + s[f]) ; 


(-1,0,1) 


transp(-l, / + 1, a; + s[f]) ; transp( — 1, /, x) ; 


end case 




end procedure. 




Figure 9: Algorithm Update_Perni. 



5 Final remarks 

Numerical evidences show that if we change the generating order of algorithm Gen_Colex 
as for Genl_Gray, but without restricting it to subexcedant sequences, then the obtained 
list for bounded compositions is still a Gray code with the closeness definition slightly 
relaxed: two consecutive compositions differ in at most four adjacent positions. Also, T. 
Walsh gives in |8] an efficient generating algorithm for a Gray code for bounded compo- 
sitions of an integer, and in particular for subexcedant sequences. In this Gray code two 
consecutive sequences differ in two positions and by 1 and —1 in these positions; but these 
positions can be arbitrarily far, and so the image of this Gray code through the bijection 
i/j defined by relation ([T]) in Section [3] does not give a Gray code for permutations with a 
fixed index. 
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